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ABSTRACT 

Pumping routes of masers can in principle be recovered from a small matrix of master 
equations, at an advanced stage of elimination, by tracing back the coefficients to a 
set of unmodified all-process rate coefficients, drawn from those which appeared in the 
original set of master equations, prior to any elimination operations. The traceback 
is achieved by logging the operations carried out on each coefficient. There is no 
guarantee that a pumping scheme can be represented as a small set of important routes 
in this way. In the present work, the traceback method is applied to a model which is 
typical of a large volume of parameter space which produces very strong inversions in 
the main lines of the rotational ground state of OH, at 1665- and 1667-MHz. For both 
lines, the pumping scheme is largely restricted to the 2 n 3 / 2 stack of rotational levels, 
and it is possible to list a comparatively small set of routes (less than ten) which 
provide more than 80 per cent of the inversion. In both cases, the strongest, and 
simplest, route consists of a radiative upward stage, to the 2 n 3 / 2 , J = 5/2 rotational 
level, followed by a collisional dc-cxcitation to the rotational ground state. 
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1 INTRODUCTION 

Sophisticated, many-level, non-local-thermodynamic- 
equilibrium (NLTE) computations are now routinely used 
to generate observables, such as the emergent flux and po- 
larization, from maser regions. Such observables are based 
on NLTE molecular energy-level populations and associated 
radiation fields, which are the fundamental outputs from 
the calculations. In general, however, computations of this 
type yield far less detailed information about the pumping 
schemes that produce the maser inversions because the 
NLTE solutions do not explicitly store pumping routes. 
The combination of radiative transfer and kinetic master 
equations that comprises a model of a maser environment 
typically loses track of which pieces of molecular data, 
from the (typically) thousands that form the input to the 
computation, are responsible for the resulting inversions. 

Our knowledge of maser pumping schemes remains 
quite patchy. In the case of OH, the pumping scheme for the 
1612-MHz line in long-period variable and supergiant stars is 
quite well understood (Elitzur, Goldreich & Scoville (|1976p : 
Elitzur (|1981[> ). Absorption of 35- and 53-/im radiation is 
required to lift population from the ground rotational state 
to the J = 3/2 and J = 5/2 rotational levels of the 2 Hi/2 
stack, followed by a series of radiative decays to the upper 
state of the 1612-MHz line. Dickinson (1987) predicted from 



IRAS data that the two pumping lines should make similar 
contributions to the pump. A discussion of modern searches 
for the pumping lines, and an application of the technique 
used in the present work (to an OH/IR star envelope that is 
far from typical) appears in Gray, Howe & Lewis (12005p . In- 
vestigations into the pumping of the OH ground-state main 
lines in stellar envelopes have also been carried out (Collison 
& Nedoluha ([1555) 1. 1(1555)1 V 

A much less complete picture emerges for OH in star- 
forming regions, particularly for the ground-state main lines 
at 1665 and 1667 MHz. One ingredient that is widely be- 
lieved to be important in the pumping of these lines is far- 
infrared (FIR) line overlap (Litvak ((1969)) : Lucas |Tl980)l : 
Bujarrabal et al. (|1980[) ). Detailed models produced after 
corrected collisional data became available (Dixon, Field & 
Zare (|1985f) ; Andresen, Hausler & Liilf i|1984[l ) have added 
further insights. Piehler & Kegel (| 1989 [) reject a photodis- 
sociative pump, whilst Kylafis & Norman (119901 1 reject a 
predominantly collisional scheme. Detailed many-level mod- 
els, for example, Cesaroni & Walmsley (11991t , Gray, Doel & 
Field (|1991|> , and Pavlakis & Kylafis (I1996j> suggest a combi- 
nation of FIR line overlap and an FIR continuum radiation 
field are significant components of the pump, but the details 
remain obscure. Recent work on the pumping of 1667-MHz 
masers in megamaser sources (|Yu 200511 links the inversion 
with FIR radiation at 60 /im. 
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An alternative approach to the analysis of pump- 
ing schemes is required: one which reveals significantly 
more detail than broad inferences from NLTE computa- 
tions. One such alternative was outlined by Sobolev (1986) 
and references thererin. This method involves treating the 
population flow as a set of cycles, with varying num- 
bers of links. These cycles were later compared by anal- 
ogy with electronic circuits governed by Kirchoff's Laws 
( |Sobolev fc Deg uchi 1994a]). Some subset of these cycles will 
be important for sustaining each maser inversion. Find- 
ing the strongest component of the pump was introduced 
as an analytic optimisation problem, but the suggestion 
that was followed is that the whole method be developed 
as a Monte-Carlo computer code, providing mainly statis- 
tical information about pumping schemes. This program 
was used to analyse the pumping scheme for water masers 
(Sobolev 1989), where it was estimated that a very large 
number of cycles would need to be traced in order to obtain 
75 per cent of the population flow. The most sophisticated 
use of this Monte-Carlo scheme is the analysis of pump- 
ing routes in methanol masers ( |Sobolev fc Deguchi 1994a] 
|Sobolev fc Deguchi 1994b] ). This work was successful in 
identifying 'bottleneck' transitions, following modifications 
due to saturation, and in computing the percentage of 
the maser flux produced as a function of the number of 
links in a cycle. A significant proportion of the flux de- 
pended on complicated cycles involving at least 25 links 
dSobolev fc Deguchi 1994a] >. 

The method I use to analyse pumping routes in the 
present work is related to the flow cycles discussed above 
( |Sobolev fc Deg uchi 1994a), but here the relationship be- 
tween simple analytic expressions for the inversion, and rate 
coefficients, from a partially eliminated matrix of master 
equations, is made explicit. The method also draws heavily 
on the simple interpretation of the decomposition of a rate 
coefficient, after an arbitrary number of matrix eliminations, 
into two sub-expressions. One of these is the same coefficient 
at an earlier stage of elimination, and the other is a route via 
a level equal to the row and column just eliminated. Further 
details of the method are explained in Gray, Howe & Lewis, 
(2005). A summary is given in Section 2 below. 



2 A SUMMARY OF THE TRACEBACK 
METHOD 

I begin with some definitions. The symbol k^, y is an all- 
process rate coefficient, representing the flow of population 
from energy level x to energy level y in a system of master 
equations which is at elimination stage p. For historical rea- 
sons, the elimination stage is defined to be p — n+ 1 where n 
is the number of rows (or columns) remaining in the matrix. 
As elimination proceeds, n decreases, and has the value 3 
when the matrix has been reduced to 2 x 2 form. I define 
the original matrix to be the matrix formed from the set of 
master equations prior to any elimination operations being 
carried out; it is a square matrix of size N x N, where N is 
the number of energy levels in the molecular model. If x = y, 
the coefficient is diagonal, and takes on the meaning of the 
sum of all rate coefficients out of level x. It can be proved 
by induction (see Appendix A) that, provided this definition 
of a diagonal coefficient holds at p = N + 1, it holds for all 



smaller values of n until matrix elimination reaches a 2 x 2 
form. Where it is necessary to raise a coefficient to a power, 
it is always enclosed in brackets with the elimination stage 
inside, and the power outside the brackets. 

The coupled radiative transfer and kinetic master 
equations are solved by a standard numerical method 
l| Jones et al. 1994f) . The same equations are then re-solved 
by a naive matrix elimination technique, in which an all- 
process rate coefficient, at elimination stage p — 1, can be 
written as 

£- p_1 — k p + k p k p lk v (11 
H,j — 31 K i,p-i K p-i,j/ K p-i,p-i W 

where the —sign applies only to modification of diagonal co- 
efficients, and where the denominator in the second term on 
the right-hand side, a diagonal coefficient, acts as a normal- 
ising factor, and this whole term represents the introduction 
of a new link to the coefficient, transfering population from 
level i to level j via level p — 1. 

Operations of the type shown in eq.{TJ are logged and 
divided into three types, depending on the relative sizes of 
the two terms on the right-hand side. If the first term is 
larger than the second by a factor of at least 1/e, where e is 
a parameter less than 1, the operation is flagged as 'unmod- 
ified'. This instructs the tracer code, introduced in Gray 
et al. (2005), to treat k v ~ x as unchanged from its previous 
value. If the first term is smaller than the second by a factor 
smaller than e, the operation is flagged as 'replacement': in 
this case the fcf . is discarded in favour of the new route via 
level p— 1. The third option (an 'amendment') requires that 
both parts of the expression be kept. 

The computer code tracer takes coefficients from a 
matrix reduced to a state where p is small (typically 4 or 
5), and expands it back to forms with higher p, using the 
information contained in the operation log to simplify the 
expressions as far as possible. Values for inversions calcu- 
lated by the standard numerical method are compared with 
those computed from coefficients returned by tracer as a 
check that e has been made small enough (or conversely 
that sufficient information has been retained). 

Providing that the pumping scheme is sufficiently sim- 
ple (and there is no guarantee of this) it is possible to use 
tracer in several stages to expand coefficients at elimina- 
tion stage p ~ 4 back to the original coefficients of the un- 
modified matrix, that is where p = N + 1. Once this has 
been achieved, the coefficient k\ 3 -, for small p, has been rep- 
resented in terms of original all-process rate coefficients that 
can be expressed directly in terms of the molecular parame- 
ters supplied to the model (Einstein A- values and collisional 
rate coefficients) and the radiation energy density, or mean 
intensity. 

2.1 inversions 

Expansion of a single coefficient does not, of course, reveal 
the pumping scheme. In order to do this, we need to re- 
cover the net effect of a set of coefficients which together 
yield an inversion in the transition of interest. The method 
used is simply to find an analytic expression for the required 
inversion in terms of coefficients at small p and then use 
tracer on all antagonistic pairs of coefficients which repre- 
sent forward (pumping) and reverse (anti-pumping) routes. 
It is important to note that some terms, which are large 
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in the tracer expansion of an individual coefficient, may 
contribute almost nothing to an inversion because they are 
paired with a term of the same magnitude in the expansion 
of the reverse coefficient. I give below the analytic formu- 
lae for the inversions in the main-line ground state maser 
transitions of OH. The formula for 1665 MHz (level 3 to 
level 1) is given in terms of coefficients with p — 4, whilst 
the analogous formula for 1667 MHz is written with p — 5. 
For a listing of level numbers in terms of the more usual 
quantum-mechanical designations for OH, see Table [T] The 
level numbers used in the present work are ordered upwards 
by energy, with level 1 being the ground state. 
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where X = fc? ifcf 3 — fci 3^3 1 and M is the total number 
density of OH used in the model. 

Equations [5] and [3] have a number of important fea- 
tures: both have, multiplying the outer bracket, a positive 
definite term which modifies the overall magnitude of the 
inversion, but which does not decide its sign. It is obvious 
that X is positive definite because the diagonal term fcf 1 
must contain k\ 3 as part of its sum, and similarly fc| 3 con- 
tains fcfi. The negative part of the expression, as written 
above, is therefore cancelled exactly. The denominator D is 
shown to be positive definite in Appendix B. Both inver- 
sions are quoted per magnetic sublevel, giving rise to the 3 
in the denominator of eq.© amd the 5, in eq.Q. Inside the 
outer bracket, groups of coefficients are written in antago- 
nistic pairs, which represent the flow of population between 
the pair of levels that form the transition of interest. The 
pair of coefficients involving only the upper and lower levels 
of the transition of interest will be referred to as forming 
the 'direct' route (even though, on expansion, many other 
levels will in general be involved) . Other pairs, which involve 
one or two other un-eliminated levels, will be referred to as 
forming 'indirect' routes, because they involve levels other 
than those that make up the transition in question before 
any expansion has been performed with tracer. The reason 
for the greater complexity of eq.([3} is that it is evaluated for 
p = 5 (a 4 x 4 matrix) whilst eq.© uses p = 4. 



3 THE OH PUMPING MODEL 

The pump-route traces for both main lines were drawn 
from the same model. This was a single computation 
chosen from a large parameter-space search using a 
slab-geometry accelerated lambda iteration (ALI) code 
(Scharmcr & Carlsson 1985). This code incorporates far- 
infrared (FIR) line overlap (Jones et al. (fl994")l : Stiff (fT992|) V 
It has previously been applied in the study of several OH 
and H2O maser environments ([Gray, Howe fe Lewis 2005[ 
|Gray 2001||Yates Field fe Gray 1997| IRandell et al. 1995|) . 



Table 1. Quantum-mechanical designations of the first twenty 
hyperfine levels of OH in Hund's case (a) notation. 



Level Number 


Designation 


1 


2 n 3/2 , j = 3/2,(1-) 


2 


2 n 3/2 ,J = 3/2,(2-) 


3 


2 n 3/2 ,j = 3/2,(1+) 


1 


2 n 3/2 ,J = 3/2,(2+) 


5 


2 n 3/2 ,J = 5/2,(2+) 


6 


2 n 3/2 , J = 5/2, (3+) 


7 


2 n 3/2 , J = 5/2, (2-) 


8 


2 n 3/2 , J = 5/2, (3-) 


9 


2 n 1/2 , j = 1/2,(0+) 


10 


2 n 1/2 , j = 1/2,(1+) 


11 


2 n 1/2 , j = 1/2, (0-) 


12 


2 n 1/2 ,j = 1/2,(1-) 


13 


2 n 1/2 ,j = 3/2,(1-) 


14 


2 n 1/2 ,J = 3/2,(2-) 


15 


2 n 1/2 ,j = 3/2,(1+) 


16 


2 n 1/2 ,J = 3/2,(2+) 


17 


2n 3/2 ,J = 7/2,(4-) 


18 


2 n 3/2 ,J = 7/2,(3-) 


19 


2 n 3/2 ,J = 7/2,(4+) 


20 


2 n 3/2 ,J = 7/2,(3+) 


Table 2. Main parameters of the ALI model 


Parameter 


Value 



Total depth 


3.0 X 10 13 


111 


Depth of chosen slab 


2.0 X 10 13 


111 


Thickness of chosen slab 


4.7 x 10 12 


111 


H 2 number density 


1.0 x 10 7 cm 


-3 


Fractional abundance of OH 


2.0 x 10 


-7 


Kinetic temperature 


30 


K 


Dust temperature 


70 


K 


Bulk velocity shift 


0.0 km s 


-1 


Microturbulence 


0.0 km s 


-1 



The only criterion for selection was that the chosen com- 
putation came from near the peak of the parameter-space 
search for unsaturated integrated gain in both the ground- 
state main lines. Fig. [T] shows unsaturated integrated gain 
plotted as a function of kinetic and dust temperatures for 
part of the parameter space search, including the model cho- 
sen. Parameters of the selected model itself are shown in Ta- 
ble(2] Unsaturated integrated gains predicted for the selected 
model were 20.039 at 1665 MHz and 35.661 at 1667 MHz. 
Saturation would limit any real maser to values in the range 
10-14 as the overall amplification factor is the exponential 
of the integrated gain through the model. Competitive gain 
in saturation ( |Field fe Gray 19 88), is effective at convert- 
ing initial inversion at 1667 MHz to maser amplification at 
1665 MHz, so the model is probably consistent with obser- 
vations in Galactic star-forming regions which show that 
1665 MHz is usually the dominant main-line, and ground- 
state, OH maser (for example, Gaume & Mutel (1987)). 
An additional inversion was present in the ground state at 
1720 MHz, but the only significant excited-state inversions 
in this model were the 2 n 3 / 2 , J = 5/2 main lines at 6035 
and 6030 MHz. 

The tracer method can only work under one set of con- 
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200 1665 MHz 




Figure 1. Integrated gains in the 1665 (top) and 1667-MHz 
lines for a dataset, which includes the chosen ALI model, plotted 
against the kinetic and dust temperatures. Physical conditions 
not shown on the axes are those from Table [2] Saturation of the 
maser transitions is not included. 



ditions at once, including the radiation field. Therefore, the 
naive matrix elimination method was run on a single se- 
lected slab drawn from the total of 85 in the model. These 
are arranged logarithmically, such that the depth of slab k 
is given by 

z k = z\{zm/zi) k (4) 



where is the depth of layer k, and there are M layers 
altogether. The selected slab was k = 80; this slab was 
chosen to have the peak inversion found anywhere in the 
model at 1665 MHz. It also turns out that the maximum 
1667-MHz inversion was also found in the same slab. The 
absolute inversions in the two lines for the chosen slab are 
A/53,1 = 7.311 x 10" 3 cm" 3 and A/9 4 ,2 = 8.683 x 10" 3 cm" 3 . 
The overall number density of OH, from Table(2] is 2.0 cm" 3 . 
The tracer analysis which follows, if it is to have any gen- 
erality, therefore includes the assumptions that the pumping 
schemes in other slabs are generally similar to those analysed 
here, and that neighbouring models, such as those forming 
the grids in Fig. [1] also have pumping schemes broadly simi- 
lar to that of the chosen model and slab. This assumption is 
probably reasonable, given that the model slab is uniform, 
and that the integrated gains in Fig. [T] vary smoothly in the 
vicinity of the chosen model. 

3.1 molecular data 

The molecular data supplied as input to the ALI code 
comprised energies of the 48 hyperfine levels used and 
Einstein A-values for the radiatively allowed transitions 
l|Destombes et al. 1977[) . complemented by collisional rate 
coefficients from Offer, van Hemmert & van Dishoeck ( 1994). 
This set of coefficients allows for collisions of OH with both 
ortho- and para-hydrogen. Molecular hydrogen was assumed 
to be distributed between the ortho- and para- species at a 
ratio of 3 to 1. 



4 THE 1665 MHZ PUMP 

As the highest level involved in the 1665-MHz transition is 
level 3, it is convenient to begin from a 3 x 3 matrix, which 
yields the inversion in eq.©. For this transition, eq.© has 
two antagonistic groups of rate coefficients which control the 
inversion: the first is the 'direct' route, ft* 3 — ^3 i an( ^ ^ ne sec " 
ond, the 'indirect' route via level 2, (^^2^2,3 — &3,2&2,i)/^2,2- 
On examining the magnitudes of these two expressions, the 
direct route was found to have the value 1.05 x 10" 4 s" 1 , 
and the indirect route, 9.07 x 10~ 5 s" 1 . Both terms are pos- 
itive, and therefore inverting, and similar in magnitude, so 
they provide roughly comparable contributions to the over- 
all pump. It is therefore clear that both routes must be taken 
into account when tracing back towards earlier elimination 
stages. Fortunately, the trace remains resonably simple, with 
two components of the direct route supplying 83% of it, and 
three components of the indirect route contributing 80.5% of 
the total in that system. This total of five dominant terms, 
which together provide 81.8 per cent of the 1665 MHz inver- 
sion, are discussed below in detail. 

4.1 Route 1 

Route 1 is the strongest term in the expansion of the direct 
pump, kf 3 — k$i. It has a rate of 6.01 x 10~ 5 s" 1 (rela- 
tive strength 1.0) and the expression in rate-coefficients at 
elimination stage 6 is (ft? sfcf 3 — ^3,5^5,1)7^5,5- Route 1 is 
the strongest of the five dominant routes, and also the sim- 
plest, because further expansion of the coefficients in the nu- 
merator with tracer produces little additional complexity. 
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Writing unmodified coefficients (at elimination stage jV + 1) 
without superscript, I found kf 5 = fei,5 and similarly for 
fc|,i- Th e expansion of fcf 3 and fcf^ produced original coef- 
ficients, and extra terms (via a list of 'amendment' opera- 
tions) . However, all but one of these were found to be weak 
and anti-inverting. The remaining route was a weakly invert- 
ing pump via level 15. However, it was negligible compared 
to the route using the original coefficients, and to Route 3 
(see below). Denoting Route 1 by Rl, the pump route can 
therefore be expressed as 



(8) 



-Rl = (fcl, 5^5,3 — ^3,5^5,1)7^5,! 



(5) 



Note that the denominator in eq.([5} is left evaluated at elim- 
ination stage 6 for simplicity. As a sum of rate-coefficients, 
it is positive definite and cannot control the sign of the ex- 
pression. Important original, or unmodified, all-process rate- 
coefficients appear in Table [3] Route 1 is also displayed by 
the black arrows in Figure [2] 

4.2 Route 2 

Route 2 is the strongest term in the expansion of the indirect 
pump. It has the expression, 

R2 = (fci j5 fcf j2 &2,4&4,3 — ^3,4^4,2^2,5^5,1)7(^2, 2^4,4^5, 5) (6) 

and a rate equal to 4.86 x 10 -5 s _1 . It has a relative strength 
(compared to Route 1) of 0.808. Unlike Route 1, the expan- 
sion of Route 2 with tracer introduces a complicated web 
of routes, not all of which have been traced fully. Three 
terms control 73 per cent of Route 2. When fully expanded 
back to unmodified coefficients, these three components of 
Route 2 are 

fcl, 5^5, 2^2, 6&6, 4&4, 7^7, 3 — fc3,7fc7 ,4^,6^^2,5^,1 



R2= 



+ 



+ 



£.4 £.5 £.6 £,5 £.8 

""1, 2^4, 4^5, 5 ""6, 6^7, 7 

fcl,5fc5,2fc2,6fc6,4fc4,8fc8,20fc20,7fc7,3 ~ ^2_B 
h A h 5 Zr-6 P P £-9 £-21 

fcl,5fc5,2fc2,10fcl0,14fcl4,4fc4,7fc7,3 ~ ^2C 
£,4 £.5 £.6 £.5 £.15 £.8 
' t 2,2' t 4,4' t 5,5' t 10,10' t 14,14' i 7,7 



(7) 



where fi2B and SI2C are reverse routes, formed from a prod- 
uct of rate coefficients with each pair of levels reversed with 
respect to those of the forward route; Q.2A is written out 
in full. The lines of eq. (JTJl correspond to the A, B and C 
subroutes in Fig. [2] The loss of 27 per cent of Route 2, miss- 
ing from eq.([7|, arises from ignoring k\ 6^6,2/^6 6 m favour 
of £5,2 in the expansion of k%2 (and a similar approxima- 
tion in the reverse route). This omission is comparable, in 
magnitude, to ignoring the fourth and fifth dominant terms 
(the remaining two from the indirect pump). The second 
and third lines in eq.© are of comparable strength, and 
these are both about 1/5 of the strength of the route on 
the first line. The web of transitions comprising Route 2 is 
represented by blue arrows in Fig. [2] (except for transitions 
already marked as part of Route 1). 

4.3 Route 3 

Route 3 comes from the direct part of the pump, and has a 
rate of 2.72 x 10~ 5 s _1 (relative strength 0.453), placing it 
third in strength of the five dominant terms. The expression 
for _R3 is 



The expansion of Route 3 leads to a complex web of routes 
most of which, for brevity, have been omitted here. The 
strongest three routes, when fully traced back to unmodified 
coefficients, appear in eq.@. These three together account 
for only 50 per cent of Route 3. Six terms appear in eq.® 
because some of the expansion involves amendment opera- 
tions which add terms. The first four terms correspond to 
Route 3A, and the last pair are Route 3B and Route 3C, 
respectively. Reverse routes are written in full except for in 
the third line. 



R3-- 



+ 



+ 



+ 



fcl,9fe9,4fe4,7fc7,3 — kzjk7 ,4,ki,$kg,l 

P P £.10 
f M,4 fi '7,7 ft '9,9 

fcl,9fcg,4fc4,8fcs,20fc20,7fc7,3 — fc3,7fc7,2ofc20,8fc8,4fc4,9fc9,l 

P P P £.10 £.21 
"•4,4"'7,7"'8,8' t 9,9' t 20,20 

fcl,9fc9,13 fcl3,10fcl0,14fcl4,4fc4,7fc7,3 — Q'iAiii 

P P £.10 £.11 p4 p5 

"'4,4 ft '7,7' 1 '9,9' 1 '10,10' i 13,13'n4,14 

fcl,9fc9,13fcl3,4fc4,7fc7,3 ~ fc3,7fc7,4fc4,13fcl3,9fc9,l 

P P £-10 £.14 
fi, 4,4 fi '7,7' i '9,9' i 13,13 

fcl,10fcl0,14fcl4,4fc4,7fc7,3 ~ fc3,7fc7,4fc4,14fcl4,ipfcl0,l 
£.5 £.8 £.11 £. 
"'4,4"'7,7"'10,10' 1 ' 

ftl,5fc5,4fe4 > 3 — fc3,4fc4,5fc5,l 



15 

14.14 



£.5 £.6 

re 4,4 re S,6 



(9) 



4.4 Route 4 

Route 4 comes from the indirect part of the pump; it has a 
rate equal to 

R4 = (fcl,5fc5,2fc2,6fce,3 - fc3,efc6,2fc2,5fc5,l)/(fc2,2fc5,5fc6,6) (10) 

with a numerical value of 1.31 x 10 -5 s _1 (relative strength 
= 0.218). 



4.5 Route 5 

Route 5 is also part of the indirect pump; it has a rate equal 
to 



R5 = 



fcl,5fc5,2fc2,lofcio,3 — fc"lofclO,2fc2,5fc5,l 

£.4 £.6 £.11 
fi '2,2 fi '5,5' i '10,10 



(11) 



with a numerical value of 1.13 x 10 5 s 1 (relative strength 
= 0.188). 

4.6 Summary 

The 1665 MHz pump, though complex, is simple enough to 
allow a large percentage of its strength to be represented in 
a modest number of terms. The overall degree of complexity 
is similar to that found for the 1612 MHz line in the OH/IR 
stellar envelopes studies in Gray, Howe & Lewis (2005}. In 
contrast to that case, in which transfer to the 2 ^li/2 stack 
from the ground state by 53 /im radiation is amongst the 
most important processes, the 1665 MHz pump tends to stay 
largely within the 2 n 3 / 2 stack, with transfer to 2 Tli/2 ap- 
pearing only in weaker terms of Route 2, and in Route 3. 
The strongest pump route of all, Route 1, is extremely sim- 
ple, comprising only two parts: a radiative FIR absorption 
from the ground state up to 2 n 3/ <2, J = 5/2, and a collisional 
decay back to level 3 in the upper part of the ground-state 
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Figure 2. Principal pumping routes for the 1665 MHz maser 
under the conditions of Table [2] Energy levels are numbered fol- 
lowing the scheme in Table[T] and are marked with a parity assign- 
ment (+ or -). Energy gaps between levels are not drawn to scale. 
Route 1 transitions are shown in black; additional transitions ap- 
pearing in Route 2 are drawn in blue, and further additions from 
Route 3 are in red. Transitions appearing from an 'A' sub-route 
have solid lines, those from a 'B' sub-route are marked by dashed 
lines, those from a 'C subroute, by dotted lines. For simplicity, 
only the forward (pumping) routes are shown for each antago- 
nistic pair, and transitions appearing in more than one route are 
drawn only once. A solid arrowhead indicates a transition allowed 
for electric dipole radiation; a hollow arrowhead indicates a ra- 
diatively forbidden transition in which population transfer can be 
considered to proceed via collisions only. 



lambda doublet. This downward leg must be collisional be- 
cause the upper part of the ground-state lambda doublet and 
the lower part of the excited-state lambda doublet have the 
same parity, making this transition radiatively forbidden. 

Most of the transitions shown in Fig.[2]are radiative, but 
there are a significant number of important links which must 
be predominantly collisional, either because the transitions 
involved are within lambda doublets and have a very small 
A- value, or they are radiatively forbidden by parity or other 
quantum- mechanical selection rules. No part of the pumping 
scheme involves levels higher than those in the 2 n 3/ / 2 , J = 



7/2 lambda doublet. 



5 THE 1667 MHZ PUMP 

This pumping scheme was traced from a 4 x 4 matrix, with 
the inversion given by eq.©. The inversion expression is 
therefore considerably more complicated than in the case 
of 1665 MHz, but the compensation for this is reduced com- 
plexity in the traceback. Seven contributions to the inversion 



were found which had an inverting effect of at least ten per 
cent of the strongest. Of these seven routes, the strongest 
pair account for 56.5 per cent of the total inversion. In fully 
traced-back form, the seven routes become: 



Rl 



^2,6^6,4 — ^4,6^6,2 



k 7 
^6,6 



(12) 



where Rl is taken to have a relative strength of 1.0, and 

fe,5fc5,lfcl, 5^5,3^3,7^7,4 ~ ^4,7^7,3^3,5^5,1^1, 5^5,2 



R2 



[^5,5] 2 ^7,7^ 
fe,5fc5,lfcl, 5^5,3^3,7^7,20^20,8^8,4 



Sir 



[£.6 121.8 £9 1.21 v 
I 5 5J 7 7 8 8 ? n *>ri 



R3 = 



2,10^10,14^14,4 



20,20 
ki 14^14,10^10, 2 



£.11 1,15 
"'lO, in"- 



0.413i?l(13) 
= 0.290fll (14) 



10 "'14,14 



R4 



&l,l(fe, 6^6,3^3,7^7, 4 — ^4, 7^7,3^3, 6^6, 2 
^7,7^8,8^ 

k x 1(^:2,6^6,3^3,7^7,20^20,8^8,4 



^^[^f.sl^io^o^ 



0.255i?l (15) 



R5 = 



+ 



fcl, 1(^2,10^10,3^3, 7&7, 4 — ^4, 7^7, 3^3, 10^10,2) 
k 8 k U X 

" , 7,7'no,io yl 

k\, 1(^2,10^10,14^14, 3^3, 7&7, 4 — fifls) 

P pi p5 V 

" , 7,7' 1 '10,10 ft '14,14 yl 



= 0.228i?l (16) 



R6 - 



R7 



fcl. 3 (fc2, 5^5, 1^1,5^5, 4 — ^4,5^5, 1^1,5/15,2) 



0.203.R1 (17) 



+ 



+ 



&3, 3(^2,5^5, 1^1,9^9, 4 — ^4,9/09,1^1, 5/25,2) 

k G k 10 X 

"■5, 5^9, 9 A 

k-l, 3 (&2, 5&5, 1^1,9^9, 13&13, 10^10,14^14, 4 — Qr7b) 

P i-10 £.11 p4 £15 V 

fi '5,5 ft '9,9 ft '10,10"'13,13'''14,14 yl - 

&3, 3(^2, 5/C5,lfcl, 9/C9,13fcl3, 4 — QR7C) 



£.6 £.10 £14 Y 
ft 5, 5^9, 9^13, 13^ 



= 0.112i?l 



(18) 



where some reverse routes have been compressed to the 
Sl-notation, as for 1665 MHz. For the equations in the set 
eq. (|13[l - eq. (|15[l where there is more than one line, the lines 
are ordered by inverting strength. The first line is termed the 
'A' route, then the 'B' route. Only Route 7 has a 'C route. 
In all these equations, inverting strengths have been given 
relative to the expression in eq. (|12[l . The absolute value of 
this expression for Rl is 6.89 x 10 _5 s _1 . The seven routes 
represented by eq. (|12[) - eo. (|15fl are shown on a schematic 
energy-level diagram of OH in Fig. [3] 



5.1 Summary 

The 1667- and 1665-MHz pumps share many features: the 
strongest pair of routes are confined to the 2 n 3 / 2 stack of 
levels, the strongest route does not involve levels higher than 
the 2 n 3 / 2 , J = 5/2 rotational state and the weaker routes 
make much more extensive use of the 2 Ui/2 levels. The 
similarity is particularly striking for the strongest pumping 
route in each line: for both lines, this consists (consider- 
ing the forward route only) of a radiative absorption from 
2 n 3 / 2 , J = 3/2 to 2 n 3 / 2 , J = 5/2, followed by a collisional 
de-excitation back to 2 H^/2,J = 3/2. For 1665MHz, the 
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Table 3. Radiative and collisional rate coefficients for 1665 MHz 



Figure 3. Principal pumping routes for the 1667 MHz maser 
under the conditions of Table[2] Markings as for Fig.[2]except that 
the colour code is extended to seven colours: transitions appearing 
in Route 4 are drawn in green and then Route 5 in gold, Route 6 
in pink and Route 7 in turquoise. 



route, in levels, is 1 — > 5, followed by 5 — > 3. The analo- 
gous route at 1667 MHz is 2 -> 6, followed by 6 -> 4. A 
similar analogy can be drawn for the second-strongest route 
in both lines, where a collisional de-excitation is bracketed 
by a pair and a triplet of radiative lines. The first group 
operates between the lower halves of the 2 Tl 3 / 2 ,J = 3/2 
and 2 n 3 / 2 , J = 5/2 lambda doublets, and the second group 
links the upper halves. The swap from the lower to the upper 
halves is allowed by the intervening collisional de-excitation. 



6 THE PUMPS AT A DEEPER LEVEL 

So far, the present work has derived pumping routes for OH 
masers in terms of all-process rate-coefficients. Therefore, 
we know what routes are responsible for most of the steady- 
state inversion, but some details of the underlying physics 
are still missing. To proceed, we can expand the all-process 
rate coefficients presented in eq.© - eq. (|ll[) and eq. (ll2|l - 
eq. (|18[) in terms of their radiative and collisional parts, so 
that the importance of routes can be explained in terms of 
radiation fields and molecular parameters. Some of the ex- 
pressions are very complex, so I present here the analysis for 
the most important route at 1665 MHz, and for its analogue 
at 1667 MHz, where the analysis of the pumping schemes 
leads to a very simple physical understanding. 



Quantity Value (Hz) 





1.24 x 10" 


i 


B 8 ,iJi,s 


2.16 x 10" 


2 




3.60 x 10- 


2 


Cs,i 


4.35 x 10" 


5 


Cl,5 


1.31 x 10" 


6 


C5,3 


3.75 x 10- 


4 


C3,5 


1.13 x 10- 


5 



6.1 1665 MHz 



The strongest component of the 1665 MHz pump is given by 
the expression in eq.([5]l. The absolute value of the inversion 
it generates requires multiplication of Rl by the constants 
outside the main brackets in eq.®- Here, the only term of 
interest is the antagonistic part of eq.([5j), so we group the 
denominator with the constants and expand the expression, 



y = fcl,5fe,3 — &3, 5^5,1 



(19) 



As the 5^3 transition is radiatively forbidden, expansion 
of eq.(|19|> in terms of the radiation field and molecular pa- 
rameters yields 

y = (Bi, 5 Ji,s + Ci, 5 )C 5 ,3 - C 3 , 5 (A 5 ,i + Ji, 5 B 5 ,i + C 5 ,i)(20) 

where A x , y , {B XtV ) are Einstein-A (B) coefficients and the 
C x ,y are first-order collisional rate coefficients for popula- 
tion transfer from level x to level y. The values of the seven 
important quantities appearing in eq. (|20p are tabulated in 
Table [3] It is easy to see from Table [3] that there is a hierar- 
chy of these rate coefficients. The stimulated emission and 
absorption terms are about five times smaller than the spon- 
taneous emission rate, but of order 100 times larger than the 
collisional rate coefficients. A sensible initial simplification 
of eo. (|20p is therefore to ignore terms that are products of 
collisional coefficients, leaving 



-Bi.sJi^Cs.s — 63,5(745,1 + Ji, 5-65,1) 



(21) 



Another glance at Table [3] shows that 6*3,5 is smaller than 
6*5,3 by a factor of about 30. Ignoring the stimulated emis- 
sion compared with the absorption reduces eo. (|21|l to 



y ~ Bi, 5 Ji, 5 C 5 ,3 - C 3 ,sA 5 ,i 



(22) 



The Einstein B-value and the upward collisional rate- 
coefficient can be expressed in terms of the downward ex- 
pressions, producing 



V = (C5,3ff5/fll)[-B 5 ,lJl 



A 5 ,ie" 



hv 



(23) 



where g\ and 35 are the statistical weights of levels 1 and 5 
respectively, v is the line-centre transition frequency of the 
5 — » 1 transition, and Tk is the kinetic temperature. Express- 
ing the Einstein B-coefficient in terms of the A-coefficient, I 
obtain 



y = 



A5,lC 5 ,3ff5 



1) 



Ji 



B V {T K ) 



l + e 



(24) 



where B v (Tk) is a Planck function at the local kinetic tem- 
perature (for the chosen slab in the ALI model). The group 
of terms outside the square braket in eq. (|24fl helps to set the 
overall speed of the pump, but does not control its direction. 
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For the moment, group these terms as K, and noting that 
the exponential in eq. (|24[) is vastly smaller than 1, the final 
form for y is 



y~K[{Ji, s /B v {T K ))-l] 



(25) 



Two things are responsible for the effectiveness of the 
pump in eq. (|25[) . The first is that the mean intensity of radi- 
ation in the 5 — > 1 transition must be greater than a Planck 
function at the same frequency, based on the kinetic tem- 
perature. The second is that the energy gap between the 
2 n3/2, J = 3/2 and 2 n 3 /2, J = 5/2 rotational levels is large 
compared to \zTk'- this condition makes the upward rate co- 
efficient in the collision-only 5 — > 3 transition much weaker 
than its downward companion. 

The first of these requirements - that the mean intensity 
exceed the black-body function at the kinetic temperature 
- can be explained in terms of three physical processes, ei- 
ther singly or in combination. The first is the presence of a 
radiation field at a temperature higher than Tk- The ALI 
model has dust at a temperature T d = 70 K, so continuum 
radiation is a possible source of the necessary mean inten- 
sity. The second possibility is that radiative transfer effects 
drive the mean intensity to a level well above what would 
be expected for LTE at 30 K, the kinetic temperature. The 
third possibility is line overlap with another FIR transition. 
However, in the model used here, the 5 — > 1 and 6^2 
transitions are not members of the overlapping groups, so 
the dust continuum and/or NLTE radiation transfer must 
be responsible. 

For the 1665 MHz pump, I compare the rate -85,1/5,1 
from Table [3] with three other relevant radiative rates. The 
first of these is the LTE rate, B 5 , 1 B V (T K ), which is 2.28 x 
10~ 3 Hz for T K = 30 K. The actual rate in the 5 -> 1 line 
exceeds this by a factor of ~ 10. The second is the black- 
body rate at the dust temperature, B$ : iB v (Td) = 2.70 x 
10 -2 Hz - substantially larger than the actual rate. Thirdly, 
there is the rate -Bs,iSi,5 generated by the source function 
in the chosen slab, which has the value 1.97 x 10 -2 Hz. The 
actual value exceeds this slightly, so the actual radiative rate 
found in the transition must depend on both a substantial 
optical depth in the 5 — > 1 line, and the presence of the 
dust continuum. This view of the inversion mechanism is 
consistent with Fig. [T] where there is little or no inversion 
for models where Td < Tk- 

It is, in fact, possible to give a reasonably detailed phys- 
ical explanation of the radiative part of the pump, that is to 
explain eq, (|25|l in terms of the details of the radiative trans- 
fer. The ultimate source of the large mean intensity in the 
5 — * 1 line is the boundary condition on the inner (high op- 
tical depth, and the more remote from the observer) bound- 
ary of the model. This specifies that the continuum becomes 
optically thick, and that the abundance of OH falls to zero. 
As a result, the radiation field at the inner boundary is a 
Planck function at the (constant) dust temperature, whilst 
the mean optical depth in the line tends to a large, but fi- 
nite, value. In the model analysed in the present work, the 
mean line optical depth reaches tm = 265. All the slabs 
which make up the numerical model (where r < tm ) will 
be referred to here as the line zone. The line zone includes 
the slab analysed with tracer. Although the line zone also 
contains dust at the same temperature as the boundary, it 
is too optically thin to contribute significant radiation to 



the pumping line. Therefore, it is definitely the boundary 
which is responsible, and any simplified model can assume 
that opacity is provided by the line only in the line zone, 
and by the continuum only in the boundary. High line opti- 
cal depths (> 150 ) in the analysed slab and those nearby, 
suggest that a diffusion approximation may lead to a rea- 
sonable physical understanding of the line zone close to the 
inner boundary. 

If a radiation diffusion approximation is assumed (see 
Appendix C), an analytical solution to the radiation transfer 
problem can be found. Within the restrictions discussed in 
the Appendix, the mean intensity in the line zone is given 
by, 



J(t) = B V (T K ) + {B v {Td) - B v (T K ))e 



K(t m -t) 



(26) 



Equation [26] has sensible limits: the mean intensity tends 
to a Planck function at the dust temperature as the opti- 
cal depth approaches the inner boundary, whilst at small 
optical depth, the limit is a Planck function at the kinetic 
temperature. Of course eq. (|26|l has no semblance of valid- 
ity for small values of r. Re-arranging eo. (|26|l to look like 
eg. ((25)l yields, 



At) 
B v (Tk 



- 1 = 



B v {T d ) 
B U (T K ) 



- 1 



-\[$C,{t u -t) 



(27) 



so in the model used in the present work, the boundary 
value of the bracket ( 10.87 ) has decayed to 8.496 over 
an optical depth shift of 32.4. This shows that the cru- 
cial parameter is £. If it were not present in eq. (|27[) , the 
mean intensity would decay to a value dictated by the ki- 
netic temperature well within the first slab of the line zone, 
and no pumping radiation would be transported far from 
the inner boundary. The full definition of £, the scattering 
parameter, is given in Appendix C for a two-level system. 
However, in the case of the 5 — > 1 line in the present work, 
with parameters given in Table [3] it is well approximated by 
£ ~ C 5>1 /A 5>1 = 3.51 x 10" 4 . This results in an effective op- 
tical depth scale thinner than that based on the line-profile 
mean by a factor of \/3C = 0.032. This scale does not match 
the decay in J found in the model, but this is not surpris- 
ing given that the diffusion approximation is not strictly 
valid, and relies on the assumption of a two-level system. 
The important point is that the effective optical depth scale 
is vastly thinner than the line scale because of a small value 
of £. The small value of £ in turn depends upon the fact 
that the Einstein A-value in the 5 — > 1 line is vastly larger 
than the downward collisional rate coefficient. Any radia- 
tion absorbed in this line has a high probability of being 
re-radiated, in the same line, rather than being collisionally 
de-excited, which would lead to thermalization of the en- 
ergy. A large fraction of the absorption therefore behaves 
effectively as scattering, so that the effects of the boundary 
continuum on pumping can be felt in the slab studied with 
tracer, and indeed substantially closer to the observer. 

It is perhaps worth noting that photons travelling di- 
rectly from the boundary do not actually penetrate very far 
into the model. This radiation also does not decay expo- 
nentially as tm — t, because the wings of the Gaussian line 
profile are always substantially optically thinner than the 
line mean. Very roughly, the actual inner boundary radia- 
tion penetrates as 1/(tm — r), for optical depth shifts > 1, 
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so it would have only a marginal effect on the slab stud- 
ied in the present work, when compared to the 'effectively 
scattered' radiation discussed above. 

An interesting rider to this investigation of the leading 
1665-MHz pump is to consider the symmetry between the 
routes linking levels 3 and 1 via level 5, as discussed above, 
and those linking levels 3 and 1 via level 7. If one assumes 
similar Einstein coefficients and collisional rate coefficients, 
one would conclude that a radiative excitation from level 3 to 
level 7, followed by a collisional decay from 7 to level 1, would 
provide an anti-inverting mirror to the similar transitions 
via level 5, negating the inverting effect of the latter. The 
reason that this does not happen in the model discussed here 
is that the downward rate-coefficient, 67,1 = 5.00 x 10 _J Hz, 
is considerably smaller (by a factor of 0.133 ) than (75,3 (see 
Table[3}- Importantly, 6*5,3 is the largest term in the tracer 
expansion of k% ,3, which follows from eq.©, and the fact 
that the 5 — >■ 3 transition is radiatively forbidden. Therefore, 
direct collisional transfer in this route exceeds all indirect 
methods of transfer between levels 5 and 3, and the large 
difference between Cs,3 and 6*3,5 is significant in creating 
inversion. By contrast, 6*7,1 is only the third most important 
route in the expansion of fcf j. It is exceeded in importance 
by radiative routes transfering population from level 7 to 
level 1 via levels in the 2 U-i/2 stack. Therefore, compared 
with 5^3, direct collisional transfer between level 7 and 
level 1 plays a much less prominent role, and the difference 
between 61,7 and 67,1, though large, is unimportant as an 
anti-inverting mechanism. 

At a still deeper level, it is possible to explain why the 
collisional rate coefficient 67,1 is so much smaller than 6*5,3: 
both coefficients are constructed from contributions due to 
collisions with ortho- and para-H2. The rate coefficients from 
collisions of OH with ortho-H2 are similar for both transi- 
tions, with the coefficient for 7 — > 1 stronger by a factor of 
1.08. However the rate coefficients for the two transitions 
derived from collisions with para-H2 are very different. This 
'parity propensity' of the OH plus para-H2 collision cross 
sections was noted by Offer et al. (|19941 l. In fact, the ratio 
of the rate coefficients for the 5 — > 3 and 7 — »■ 1 transitions, 
due to collisions with para-H2, is 8.37 at the 30 K kinetic 
temperature of the model. The strong parity propensity of 
the collisions with para-H2 is apparent at this temperature 
because the equilibrium abundance of the ortho-H2 is small 
(about 2.5 per cent). At higher temperatures, the abundance 
of ortho-H2 would be expected to rise, and the efficiency of 
the pump to fall, as the ortho- to para-H2 abundance ratio 
tends towards its high-temperature value of 3. A glance at 
Fig.[T]suggests that this is indeed the case, with the pumping 
mechanism becoming ineffective below a kinetic temperature 
of 100 K, where the ratio of the contributions of the ortho- 
and para-hydrogen species to the overall rate coefficient is 
1.5. 



6.2 1667 MHz 

The strongest pump route at 1667 MHz can be analysed in 
a very similar way to the 1665 MHz route discussed above. 
The inversion expression to be expanded is 



Table 4. Radiative and collisional rate coefficients for 1667 MHz 



Quantity Value (Hz) 



M,2 _ 


1.38 X 10" 


-1 


B%, 2^2,6 


2.31 x 10" 


-2 


B2, 6-^2, 6 


3.24 x 10" 


-2 


C*6,2 


5.16 x 10" 


-5 


C*2,6 


1.30 x 10" 


■6 


C*6,4 


4.74 x 10" 


-4 


C*4,6 


1.20 x 10" 


-5 



and the version of eq. (|28[) fully expanded in terms of molec- 
ular parameters and the radiation field in transition 2^6 
is 

y = (£2,6 J 2 , 5 + 6 2 , 6 )6 6 ,4 - 64,6 (A 6 ,2 + J2, e B 6 , 2 + 6 6 , 2 ) (29) 

Values of the various terms appearing in eq. (129[l appear in 
Table |4] Following through the same steps as for 1665 MHz, 
the final expression for y' is 



y' ~ K'[(J2,8/B v (Tk)) - 1] 



(30) 



It can be shown that the mean intensity in the 2—^6 line 
which pumps 1667 MHz is actually smaller than the corre- 
sponding mean intensity in the 1 — > 5 transition. On this 
basis y' in eq. (|30[) would be expected to be smaller than y 
for the case of identical constants K = K' . The reason that 
the 1667 MHz inversion is larger, per sublevel, by a factor 
of 1.188 therefore requires that either K > K' , or that the 
inversion ratio is explained by differences in the external 
constants appearing in eq.(j2} and eq. (J3jl . If we take the in- 
version per sublevel, then the ratio of the external constants 
is 



Ap 4 
Ap 3 



3(&i,i&3,: 



K 1,3 K 3,: 



ei.4 u5 
0^2,2^4,4 



(31) 



When evaluated, R ex t = 1.093. Therefore, part of the reason 
for the larger inversion at 1665 MHz lies in the ratio of the 
constants, R c = K' / K. Ignoring the exponential terms, the 
ratio of which is 1.0 to better than one part in 10 4 , R c is 
given by 

^6,266,4<?6<?l 



Rc = 



(32) 



y 



^2,6^6,4 — ^4,6^6 



(28) 



^45,165,3(7235 

which evaluates to R c = 1.1794. One reason why the in- 
version at 1667 MHz is larger than the corresponding inver- 
sion at 1665 MHz is therefore that both transitions in its 
strongest pumping routes are faster than their analogues 
in the 1665 MHz pump. For comparison of the Einstein A- 
values and collisional rate coefficients, see Tableland Ta- 
ble [3] The radiative part of the pump for 1667 MHz depends 
on the efficient transfer of radiation from the inner boundary 
of the model by effective scattering of radiation in the 6^2 
line in a very similar manner to the detailed discussion for 
1665 MHz, given in Section 6.1. 

The similarity of the 1667- and 1665-MHz pumps also 
extends to the symmetry breaking between the pumping 
route (level 2 to level 4 via level 6) investigated above, and 
its anti-inverting mirror, taking population from level 4 to 
level 2 via level 8. Just as in the case of the 1665-MHz pump, 
the inverting route has the direct collisional transfer as the 
leading term in the tracer expansion of its downward tran- 
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sition ( 6 — > 4 in this case), but the direct transfer is rela- 
tively unimportant in the downward transition ( 8 — ► 2 ) of 
the anti-inverting route. The weakness of the direct 8^2 
collisional transfer is based on the size of the collisional rate 
coefficient (C8,2/C6,4 = 0.18), which in turn depends on the 
overwhelming contribution of the para-H2 species, with its 
large parity propensity, to collisions at low temperatures. 



7 EFFICIENCIES 

I discuss here the efficiency of overall schemes and the effi- 
ciencies of important individual routes. 1 define the efficiency 
of a scheme or route as the ratio of the rate at which it 
produces inversion divided by the total population flow rate 
through the same scheme or route. For example, considering 
the overall 1665 MHz pumping scheme, there are the direct 
and indirect terms from eq.©. For the direct part of the 
pump, the efficiency is 



8 DISCUSSION 



K l,3 



^1,3 



(33) 



which evaluates to 4.1 per cent. Similarly, the indirect part 
of the pump is 3.34 per cent efficient. Weighting these two 
efficiencies by the contribution each makes to the inversion, 
the overall efficiency of the 1665 MHz pump is 3.75 per cent. 
In the case of the 1667 MHz pump, obtaining the overall ef- 
ficiency is a little more laborious, but still straightforward, 
and the mean of all the routes in eq.Q, weighted by contri- 
bution to the inversion gives an efficiency of 4.60 per cent. 
The third term in the brackets of eq.Q is notable in having 
the largest individual efficiency, of 5.13 per cent. Overall, the 
1667 MHz pump is the more efficient, which is unsurprising, 
since the inversion in this line is also the larger. 



7.1 individual routes 

Individual routes within the overall scheme can have sub- 
stantially higher efficiencies than those for the scheme as a 
whole. I define such internal efficiencies as those calculated 
for a route when the denominator is confined to considera- 
tion of that route only. For example, the strongest pumping 
route at 1665 MHz has the internal efficiency, 

&1,5&5,3 — ^3,5^5,1 



&1,5&5,3 + &3, 5^5,1 



(34) 



which has the numerical value of 78.2 per cent. The 
analagous route in the 1667 MHz system, operating via levels 
2, 6 and 4 is internally 77.6 per cent efficient. These figures 
drop to ~ 2.5 per cent when the routes are expressed as 
part of the expansion of eq. (|33|l . or its 1667 MHz analogue. 
Routes with many links are not noticeably less efficient. For 
example, Route 2B in the 1665 MHz scheme, which climbs 
to level 20 in the 2 n 3 / 2 ,>/ = 7/2 rotational level, has the 
internal efficiency 



£R2B 



fel, 5^5, 2^2, 6^6, 4&4, 8&8, 20^20, 7&7, 3 — ^2S 
fcl, 5&5, 2^2,6^6, 4&4, 8&8, 20&20, 7&7, 3 + ^2S 



(35) 



with the numerical value of 77.0 per cent. 



The results reported here and in Gray et al. {2005 ) show that 
it is possible to trace pumping schemes in OH in considerable 
detail in at least two of the common OH-maser environments 
(Galactic star- forming regions and evolved-star envelopes). 
The pumping schemes revealed so far are actually very com- 
plex, but a large proportion of their inverting strength can 
be explained in terms of a few dominant routes. An obvious 
extension of this work is to study the additional ground-state 
OH maser environments: OH megamasers and those super- 
nova remnants which support strong masing at 1720 MHz. 
Additional work needs to be done to confirm the generality 
of the schemes analysed in the present work, both in terms 
of the spatial variation within one model, and variation be- 
tween models with different input parameters. To complete 
such an analysis of many schemes in a reasonable time re- 
quires considerably improved automation of the population 
tracing process. The explanation of the radiative part of 
the pump in terms of the transport of radiation from the 
optically thick boundary to a significant part of the model 
(many slabs) by means of a modifed optical depth scale does 
at least suggest that the same pumping mechanism operates 
within a large fraction of the current geometrical model. 

The importance of parity propensity in collisions be- 
tween OH and para-H2 in the pumping scheme for both OH 
main lines means that the ortho- to para-H 2 ratio in the 
gas of the maser zone must be well biased in favour of the 
para-species for the schemes investigated here to operate ef- 
ficiently. Given that H 2 forms in the ratio of 3:1 in favour 
of the ortho-form, as dictated by the nuclear spin statistics, 
maser action via the pumping schemes discussed here would 
require time for the species ratio to move substantially to- 
wards the thermodynamic ratio via a set of proton exchange 
reactions (Flower, Pincau des Forets & Walmslcy 2006). In 
some circumstances the timescale for these reactions could 
equate to the 'switch-on' time for main-line OH masers. 

Labelling of energy levels is of course arbitrary, so ex- 
pressions similar to eq.© and eq.([3} can be formed for 
excited-state inversions without increased complexity, pro- 
viding that, say, the 2 n 3 / 2 , J = 5/2 levels were labelled 1—4 
instead of those in the 2 n 3 / 2 , J = 3/2 ground state, ft would 
then be possible to trace the common 6030 and 6035 MHz 
inversions; alternative re-labelling would allow tracer to 
be applied to any of the other OH rotational levels that sup- 
port maser action. A likely problem with the extension of the 
method to excited states is one of stability: the elimination 
process following the naive method relies on the last levels 
eliminated having large populations as a form of numerical 
pivoting. There may, however, be ways of avoiding these sta- 
bility problems by attaching the operation log to the basic 
numerical method, provided that the interpretation in terms 
of rate coefficients can be preserved. 

It may also be possible to extend this inversion-tracing 
method to molecules other than OH. In many interesting 
cases, such as the 22 GHz maser line of water, the maser 
levels lie well above the ground state, so any attempts at 
analysis would suffer from the stability problems already 
discussed in connexion with the excited rotational states of 
OH. However, there may be an even more fundamental bar- 
rier to understanding inversions in some cases. In the case 
of methanol ( |Sobolev fc Deguchi 1994a| the authors show 
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that, for the 2o — 3_i Class II maser transition of the E- 
species, 65 per cent of the maser flux results from pumping 
cycles with a number of links which exceeds 10. In these 
cases at least, it is almost certainly not possible to write 
down a manageable expression for the inversion: either the 
number of required terms would be too great, or the ex- 
pressions too complicated, or both. If the tracer method 
is to be extended to another molecule, the first one to try is 
probably SiO: it has a relatively simple level structure, and 
a maser (though not one of the strongest or most common) 
has been observed between the lowest two rovibrational lev- 
els (v = 0, J = 1 — 0). In addition, population in each vi- 
brational level tends to be concentrated towards the lower 
rotational levels, which is likely to aid numerical stability. 

Another application which is likely to prove fruitful is 
to extend the tracer method to saturating conditions. The 
semi-classical analysis by Field & Gray (|1988| automatically 
breaks the kinetic scheme into coefficients which multiply 
the maser radiation and those which do not. The former 
affect the inversion as the maser saturates, whilst the lat- 
ter generate the inversion, like those studied in the present 
work. Both sets of coefficients in models (for example Gray, 
Field & Doel •: 1!I!I2; : are generated in a very similar man- 
ner, and, for the ground-state OH masers, would introduce 
no extra stability problems. The coefficients that multiply 
the maser intensity could be studied as a function of prop- 
agation distance along the maser, showing how the increas- 
ingly rapid transfer of population across the maser levels 
introduces newly important transfer routes. 



9 CONCLUSIONS 

I conclude that it is possible to trace the population trans- 
fer routes which cause strong inversions in the 18 cm OH 
main-line masers under conditions typical of star-forming 
regions. The number of routes which need to be traced in 
order to recover the great majority of the inversion is mod- 
est - a conclusion which does not necessarily hold for the 
kinetic schemes of other molecules, or even for OH in other 
environments. 

The pumping schemes derived for the 1665- and 1667- 
MHz lines show strong similarity, and are biased quite 
strongly to the 2 n 3 / 2 stack of rotational levels, unlike the 
case of 1612 MHz OH masers in late-type stellar atmo- 
spheres. In both lines, the strongest route contains just two 
transitions, comprising, for the forward route, a radiative ab- 
sorption, lifting population to the 2 n 3 / 2 , J = 5/2 rotational 
state, followed by collisional de-excitation to the upper part 
of the ground-state lambda doublet. This inclusion of a ra- 
diatively forbidden link to switch from the lower to the up- 
per half of the lambda-doublet system is also employed by 
the next most important pumping route in both main lines. 
Anti-inverting 'mirror' transitions are not effective because 
of the parity propensity found in collisions between OH and 
para-H2. 

The effectiveness of these routes in the 2 n 3/ / 2 stack is 
based firstly upon the fact that the energy gap between the 
2 n 3 /2, J = 3/2 and 2 n 3/ / 2 , J = 5/2 rotational states is large 
compared to the energy equivalent of the kinetic tempera- 
ture, which favours collisional de-excitation over excitation, 
and secondly on the presence of a dust radiation field which 



has a characteristic temperature which is locally hotter than 
the kinetic temperature, a result of efficient transfer of ra- 
diation from the optically thick boundary by 'effective scat- 
tering'. 
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APPENDIX A: DIAGONAL COEFFICIENTS 

The aim of this section is to prove that a diagonal coefficient, 
k v - j, at elimination stage p is equal to the sum of the all- 
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process rate-coefficients from level j to all the un-eliminated 
levels, or 



k h = 



E 



(Al) 



This is true by definition for the original matrix (p = N + 1), 
but is not obviously so for any smaller value of p. 

The state of the set of kinetic master equations at elim- 
ination stage p looks like: 



iPp-i 



— 'V-l,p-2ft>-l 



^p-l,2ft>- 
%_1 iPp- 



-k p 
i,P- 

-k p 

i,P- 



-Kin 



"^-2^=0 



= 



-*? ,pi =0 (A2) 



: =0 



where is the population of level x and the conservation 
equation coefhcients, written with a star in the superscript, 
appear in the final equation. An additional equation is now 
eliminated by the naive method: the topmost equation in 
the set eq. (|A2Jl is used to eliminate p p -x, such that 



p-2 

Pp- 1 = ] ^m,p-lP"0 Ap-l,p-l 

m — 1 



(A3) 



and this population is eliminated from all the others in 
favour of the expression in eq, (lA3[) . Concentrating on the 
arbitrary equation j such that j < p — 2, the diagonal coef- 
ficient in this equation is modified to the form 



K j,j 



k p 



k 1 ' 



k p 



k v 



(A4) 

I now assume that the desired result is true at elimination 
stage p, and prove by induction that it is true at all further 
stages with p > 2, noting that at p = 3, only a trivial 2x2 
matrix remains. This assumption allows the development of 
eq.l|A4p to 



K j,j 



k v 

k v 
>- 



k h- 



1 + 



k p 

K 3,P- 



k p 



(A5) 



I now add zero to each term on the right-hand side of 
eq. (|A5|l . but in a form which allows development of the off- 
diagonal coefficients. Equation IA5I becomes. 

UP LP LP LP 

LP _j_ "i,? -1 f -1 ' 1 "j,p-l>-l,l 
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(A6) 



and the off-diagonal coefficients are now of a form which can 
be upgraded to the next elimination stage via eq. ([!]). The 



coefficient at level p combines with the positive fraction on 
each line (except the last) to form a coefficient at stage p—1, 
leaving 



LP- 1 



k p 



i- p_1 + k p 



(^p-1,1 + •■• + ^p-l,. 



+ 



k p 4- 



+ k p 



+ *£-iJ 



(A7) 



where a common factor of fc P f> _ 1 /fcp_ 1 j has been ex- 
tracted from the negative terms on the right-hand side of 
eq. (|A6|l . I now note that the bracket in eo. ()A7P is the sum 
of all the rate coefficients taking population out of level p — 1, 
and is therefore equal to the diagonal coefficient k p _ lp _ 1 , 
given the assumption about such coefficients at elimination 
stage p. The bracket therefore cancels with this coefficient in 
the denominator, leaving two coefficients of the form k p p _ 1 
with opposite signs. Cancellation of these leaves 



k p ~ 

K 3,j 



E 



LP' 1 

K 3,m 



(A8) 



which is of the form of eq. (|Al|l . but with p replaced by p— 1. 
Therefore, if the initial assumption is true for any given stage 
of elimination, it remains true for the next. As it is true by 
definition for the original set of equations, where p = TV + 1, 
then the assumption is true for all subsequent eliminations 
at least as far as p — 3. 



APPENDIX B: DENOMINATOR 

This section shows that the external denominator (D in 
eq.Q and eq.Q) is positive definite. In terms of rate coeffi- 
cients evaluated at elimination stage p = 4, the denominator 
is given by 

D = fct, 3 (fe^*lfe|, 2 + fc!,2 fc 2,*l) + ^3,1(^1,3^2,2 + fcl,2&2,3) 



+ 



kf*ikt >3 ) 



(Bl) 



where coefficients marked with an asterisk in the superscript 
(starred) come from the conservation equation, 



pi 



--N 



(B2) 



and cannot be interpreted in the same manner as the other 
rate-coefficients via eq.{T}. At p = N + 1, these starred co- 
efficients are all equal to 1.0 and subsequent actions can 
only add combinations of positive-definite rate-coefficients 
to them. They are therefore positive definite at any value of 
p. Given this result, and the proof in Appendix A, the mid- 
dle term in eq. (|Bl|l is positive definite, and will be called B. 
Using the result of Appendix A to expand fef 3 = fcf 1 + fcf 2 , 
eq. (|Bl|l becomes 



D — ^3, 1^1, 1^2, 2 + ^3,1^1,2^2,1 + ^3,2^1,1^2,2 + ^3,2^1,2^2,1 

(B3) 



+ B + ki 2 kf 3 kti 



k 4 k 4 * k 4 
"■3,2"'i,i"': 



2.3 



following the expansion of the first and third brackets of 
eq. ()Bl[l . The first, second and fourth terms of eq. (|B3[) , and 
the first term following B, are all positive definite, and may 
be combined as A. This leaves, 



Pumping of OH Main-Line Masers in Star-Forming Regions 



D = A + fc^afc^^j, +B- fcj.afcftfc^a (B4) 

but when fc| 2 is expanded in accordance with Appendix A, 
it contains £2,3, so the remaining negative term in eq. (|B4[) 
is cancelled exactly, and D is positive definite. 



APPENDIX C: RADIATION DIFFUSION 

In a zone where radiation transfer is diffusive, application 
of the Eddington approximation (see for example Rybicki & 
Lightman (1979 )) leads to a radiation diffusion equation for 
the mean intensity of the form, 

^ = 3C(J(r)-B,(r)) ; (CI) 

where B v (r) is a the Planck function at the kinetic temper- 
ature, and the source function has been elimated in favour 
of the mean intensity, using the standard expression from 
molecular kinetics, 

S(t) = £(t)B„(t) + (1 - C(r))J(r). (C2) 



I note that eq. (|C2[) implicity assumes a two- level model: in 
the many-level case, the multiplier of J cannot be repre- 
sented as 1 — £, and both this multiplier and ( depend on 
the mean intensities and molecular parameters of all other 
radiative and collisionally allowed transitions. In the two- 
level case, the parameter £ is independent of optical depth 
if the kinetic temperature is constant, and is given by 

C ^~9iCi,u (C3) 

guC u ,i — giCi, u + guA Ut i 

where u is the upper and the lower level of the transition. 
As elsewhere in the paper, C u ,i ( Ci tU ) is the downward 
(upward) collisional rate coefficient, A U: i is the Einstein A- 
coefficient and g u , gi are statisical weights. 



